Screw dynamo in a time-dependent pipe flow 
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The kinematic dynamo problem is investigated for the flow of a conducting fluid in a cylindrical, 
periodic tube with conducting walls. The methods used are an eigenvalue analysis of the steady 
regime, and the three-dimensional solution of the time-dependent induction equation. The configu- 
ration and parameters considered here are close to those of a dynamo experiment planned in Perm, 
which will use a torus-shaped channel. We find growth of an initial magnetic field by more than 3 
orders of magnitude. Marked field growth can be obtained if the braking time is less than 0.2 s and 
only one diverter is used in the channel. The structure of the seed field has a strong impact on the 
field amplification factor. Generation properties can be improved by adding ferromagnetic particles 
to the fluid in order to increase its relative permeability, but this will not be necessary for the suc- 
cess of the dynamo experiment. For higher magnetic Reynolds numbers, the nontrivial evolution of 
different magnetic modes limits the value of simple 'optimistic' and 'pessimistic' estimates. 

PACS numbers: 47.65.+a, 07.55.Db, 61.25.Mv 



I. INTRODUCTION 

The screw dynamo is a simple dynamo model that has 
been extensively studied in dynamo theory. It is based on 
an idea due to Lortz p| and Ponomarenko |J| , according 
to which magnetic field can be generated by the helical 
motion of a rigid, electrically conducting cylinder of in- 
finite length through an infinitely extended medium of 
equal conductivity. This problem implies a simple veloc- 
ity field and leads to a critical magnetic Reynolds number 
as low as Q 



77 



(1) 




where U is the (constant) longitudinal velocity of the 
cylinder, vq is its radius and 77 is the magnetic difFu- 
sivity. This dynamo model has a discontinuous veloc- 
ity profile, and will be referred to as the 'Ponomarenko 
dynamo'. More realistic models of the screw dynamo, in- 
volving continuous and hydrodynamically realistic veloc- 
ity fields were considered by several authors. References 
^] develop a very accurate asymptotic theory for the 
screw dynamo in smooth flows, which hasjaeen comple- 
mented by numerical simulations 6]. In 7], this theor 
is applied to a number of realistic flows. Reference 
extended the numerical analysis to flows fluctuating in 
time. Other time dependent screw dynamo models were 
presented in fe] in connection with the Perm dynamo ex- 
periment. Reference 10] has put the screw dynamo into a 
larger context of slow dynamo mechanisms and pd| have 
generalized the concept to non-axisymmetric flows. The 



'Electronic address: 'Wolf gang. Dobler@kis.uni-f'reiburg.de| 



FIG. f : Photograph of a water experiment showing stream- 
lines in an initially spinning torus after abrupt braking. The 
white object in the channel is the diverter, a kind of a fixed 
ship screw that makes the motion strongly helical. From 
Ref. :9i. 



nonlinear behavior of the screw dynamo for spiral Cou- 
ette flow has been investigated with asymptotic methods 
by and numerically by [l3|. 

We are interested in the screw dynamo problem in the 
context of a new experimental dynamo project The 
basic idea of this project, introduced in |l4j . is to realize 
the dynamo effect in a strongly time-dependent helical 
flow. The flow is generated in a quickly rotating toroidal 
channel after abrupt braking, and is shown, for a water 
experiment, in Fig.H This application raises new ques- 
tions concerning the screw dynamo, which have not been 
addressed in previous studies. In particular, the flow 
of the conducting fluid will be located in a closed chan- 
nel and will be supercritical during a short time interval 
only. This requires the solution of the induction equation 
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for three-dimensional geometry in a time-dependent flow. 
Moreover, in oder to understand the satm'ation of the 
magnetic field (if it occurs during the experiment), one 
has to investigate the time evolution of the fully nonlinear 
MHD equations. The experimental scheme also requires 
strong optimization of the channel (minimization of its 
mass under optimal conductivity and wall thickness). 

In this paper we investigate the screw dynamo in a 
time-dependent flow using two different methods which 
allow different questions to be addressed. First, the anal- 
ysis of the eigenvalue problem related to the case of a 
steady velocity field gives insight into the full spectral 
structure and is numerically the most efficient approach. 
Our second method — numerical solution of the three- 
dimensional, space- and time-dependent problem — is 
numerically much more demanding, but it allows us to 
investigate the full three-dimensional structure and is the 
only approach that has the perspective of tackling the 
nonlinear problem. While the final goal is to solve the 
fully nonlinear problem and to include all geometrical 
and dynamical effects, the present paper only discusses 
the linear phase of magnetic field growth in (periodic) 
cylindrical geometry. 

Most laboratory dynamo projects use liquid sodium as 
conducting fluid, which has a magnetic Prandtl number 
Pm = V /t] « 10~^. This means that, in order to achieve 
the critical magnetic Reynolds number of a few tens, one 
operates at kinematic Reynolds number of order 10^-10'', 
which is far beyond what can be numerically simulated 
even on the largest computers. Thus, we cannot solve the 
dynamical part of the problem consistently and in the 
current paper just use mean velocity profiles as inferred 
from experiments. 

The structure of the paper is as follows. Section [H] 
specifies the equations and geometry of the problem ad- 
dressed here. In Sec. IIIII we present numerical solutions 
of the induction equation for prescribed, steady velocity 
profiles similar to those found in cylindrical pipes. If the 
velocity field is axisymmetric and identical in any cross 
section through the cylinder, the problem can be reduced 
to a one-dimensional eigenvalue problem that is solved 
numerically by discretization. This approach is suffi- 
ciently efficient to allow us to scan the space of relevant 
parameters and to isolate the cases that will be most fa- 
vorable for the realization of the experiment. Section Hvl 
presents results obtained with a three-dimensional MHD 
code solving the induction problem for a velocity field 
that varies in space and time. We finally draw some con- 
clusions about the realizability of the planned dynamo 
experiment in Sec. 



II. FORMULATION OF THE PROBLEM 

The proposed Perm dynamo experiment [T^ will im- 
plement a helical, strongly time-dependent flow of liquid 
sodium (Na(i)) through a torus, surrounded by a thin 



shell of copper (Cu), the electrical conductivity of which 
is about five times higher than that of liquid sodium. The 
role of this conducting shell is to 'anchor' the magnetic 
field lines in the exterior frame, so that the shear due 
to the velocity difference between fiow and exterior can 
enhance the magnetic field. 

Mathematically, dynamo action is characterized by the 
presence of growing solutions B(x, t) of the induction 
equation 

— = V X [(v ryV In^r) X B - 77V X B] , (2) 

satisfying the solenoidality condition 

V • B = . (3) 

Here B denotes the magnetic flux density, v the velocity 
field, /io a-iid //j- are the magnetic vacuum permeability 
and the relative magnetic permeability of the medium, 
and 77 is the magnetic diffusivity of the fiuid, related to 
the electrical conductivity ct by 77 = (/io/irc)"^. The term 
77V In /i]. = Vp arises if the relative magnetic permeabil- 
ity is a function of position and gives rise to what we call 
"paramagnetic pumping" (with an effective velocity Vp) 
of magnetic flux into the regions of enhanced permeabil- 
ity. A self-consistent description would also include the 
Navier-Stokes equation 

i9v . „ s V p i X B 

— = -(v.V)v-— ^1 + :^ + 7^Av, (4) 

Ot Q Q 

together with the condition 

V • V = (5) 

for an incompressible fluid. Here p denotes pressure, g 
density, j = V x (B/z^oA'r) the electric current density and 
V is the kinematic viscosity. The equations given above 
are complemented by boundary conditions describing the 
properties of the walls (see 14]). 

In the present paper, we restrict ourselves to the kine- 
matic dynamo problem, i.e., we consider only Eqs. (O 
and using a given solenoidal velocity field v(x, t) and 
neglect the magnetic back-reaction through the Lorentz 
force. This leads to a linear problem in B and is very 
helpful to understand the evolution for weak and mod- 
erately strong magnetic fields. This approach will allow 
for optimization of the experimental apparatus in many 
respects since the kinematic growth of B over several 
orders of magnitude is necessary for the success of the 
experiment. 

The turbulent fiow of a conducting medium will give 
rise to turbulent induction effects, which can be esti- 
mated by mean-field theory. In the present paper, we 
neglect these extra terms in the induction equation Q 
and refer the reader to the discussion in Ref. [13 . 

While the curvature of the pipe can be expected to play 
a role for the geometry of the experiment (the ratio of 
outer radius R to inner radius tq of the torus being about 
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FIG. 2: Cross section of the torus/pipe. The interior of the 
pipe, r < ro, is filled with a liquid (sodium) of electrical con- 
ductivity (jfl, the solid shell tq < r < ri (made of copper) has 
a conductivity ash and is itself surrounded by an electrically 
insulating medium (i.e. (Joxt = 0). 



3), we currently neglect it by considering a cylindrical 
pipe instead of a torus. Connecting the two ends of the 
cylinder by periodic boundary conditions and setting its 
length equal to Lz — 2ttR, we obtain a reasonable first 
approximation to torus geometry. 



III. ONE-DIMENSIONAL, 
TIME-INDEPENDENT PROBLEM 

A. Equations 

Let us consider the evolution of the magnetic field 
in a given helical flow through a cylindrical, conduct- 
ing pipe with finite wall thickness d = ri—r^ (the in- 
ner radius of the pipe is tq, the outer ri), surrounded 
by an electrical insulator; Fig. |2 shows a cross sec- 
tion through the pipe. We adopt cylindrical coordinates 
(r, (f, z) and assume an axisymmetric helical velocity field 
v(r) = [Q,rLo{r),Vz{ry\. The magnetic diffusivity and 
permeability may vary as functions of radius, rj = ri{r), 
Mr = ^J■r{r)■ 

Solutions of the kinematic dynamo problem can be rep- 
resented as a superposition of individual modes with ex- 
ponential growth. Because the coefficients in the equa- 
tions depend on r only, we are looking for solutions in 
the form of propagating helical waves 



B(r, V5,z,t) = b(r)e' 



'yt-\-i{mip-\'kz) 



(6) 



where 7 appears as an eigenvalue and is in general com- 
plex. The real part of 7 determines whether B decays 
(3^7 < 0) or grows (3^7 > 0). For a given mode, there 
exists a critical value Rm* of the magnetic Reynolds num- 
ber for which 3?7 changes sign from negative to posi- 
tive. The lowest value of Rm* is the threshold for dynamo 
action. For the type of dynamo problem considered here, 
it is typically between 10 and 100 (see Fig. [S] below for 
an illustration). 



Inserting Eq. ^ into the induction equation 101, one 
gets 



jbr + i{muj + kvz)br + 



d{rbr) 
r dr 



Rm" 



■qlDbr —b 



jb^ + i{muj + kvz)b^ + 
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duj 
= r—br 
dr 



Rm" 



dr 
ri [Db^ 



2im ^ 



dr] / d{rb^) im^ 
dr \ r dr r 



where Vp — rj din / dr , and 



D 



1 d 



dr \ dr 



(8) 



(9) 



is a Laplacian-type operator. Equations ((TJ, lO are writ- 
ten in non-dimensional form: distances are measured in 
units of the radius tq, velocity is measured in units of 
the longitudinal velocity on the axis of the flow, U. The 
magnetic diffusivity r](r) is measured in units of 77fl, which 
introduces the magnetic Reynolds number 



Rm = /ioAir,flO-flroC/ 



rpU 

m 



(10) 



Consequently, in our non-dimensional units we have rj = 
1 within the fluid and = Mr.flO'fl/Mr.shO'sh in the sur- 
rounding shell. 

In Eqs. (0 and (jHJ, the longitudinal field component 
bz does not enter; once br and b^ are known, it can be 
derived from the solenoidality condition V • B = 0, which 
yields 



-r-ri^br) - —b^ 

kr dr kr 



(11) 



(note that fc ^ for all growing modes). 

The shell r^ < r < ri can be treated as a special case 
of the above. In fact, in this case cu = Vz — 0, and the 
solution of {Tj), ||SJ) is well-known and given by 



br ± ib^ = C±I„i±i{K.r) + D±K„i±i{Kr) , 



(12) 



where k = ^/k"^ + Hoaj, and C± and D± are complex 
constants determined by the boundary conditions. Here, 
Im{') and Km{') are modified Bessel functions, which 
are related to the Bessel and Hankel functions of the first 
kind by [3 



Ini{z) 



i ™J„(iz) , K„,{z) = -f' 



-'H'^\^z). (13) 



A straight-forward approach would be to use this ana- 
lytical solution and match it to the solution in the inner 
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region r < ro (requiring continuity of the tangential com- 
ponents of the electric field) and to a potential field in 
the exterior (see below). However, we decided to use a 
simpler approach and solve Eqs iQ, © in the full region 
< r < ri for the step-like profile 



r]{r) 



r < ro 

tq < r < ri 



(14) 



Here, and for the one-dimensional results presented be- 
low, we have set fi,- — 1 everywhere. 

In the insulating outer domain r > ri, the current 
density is zero, V x B = 0, and thus B can be expressed 
in terms of a scalar potential P{r, ip, z, t), 

B = -VP. (15) 

Solenoidality of B leads to the potential equation 



ld_ 

r dr 



dP 
dr 



1 d^P d^P 



= 0, 



(16) 



and due to the symmetry of the problem P(r, z, £) can 
be written in the form 



P(r, if, z, t) = p(r)e''*+*(™'^+'=^) 
Equations 116|l and (|17|l result in 
„ 1 , 



V 

Solutions of C 
form 



-V 



m 



p = 



(17) 



(18) 



that are bounded for r — > cxd have the 



(19) 



The boundary conditions at r = ri are obtained from 
the requirement of continuity of b on the outer border of 
the conducting sheU. Together with Eqs. lfTB|) . lfT7|) 
and H19|l. this leads to 



b^iri) 
bz(ri) 



\k\n Ki„ 



in) 



m Km{\k\ri) 



kri 
m 



(20) 
(21) 



where K'^^{x) = dK„i{x) / dx. Eliminating fez(ri) from 
Eqs. ((TT)l and if^ one finally gets 

hr{ri) + ribUn) - -^ + 6^(ri) . (22) 

The inner boundary conditions follow from the regularity 
of b at r = and result in 



(23) 



KiO) = 6^(0) = 0, for |m| = 1, 
6,(0) = 6^(0) = 0, for \m\ ^ 1. 

The system together with the boundary con- 

ditions (pn| . (PI and is a non-selfadjoint eigenvalue 
problem. Dynamo action implies the existence of eigen- 
values 7 with 3^7 > 0. To obtain numerically the eigen- 
values and eigenfunctions, we replace the derivatives by 
their finite-difference counterparts, using 200-800 grid 
points for this discretization. The resulting matrix eigen- 
value problem is solved using the QR-algorithm. 



B. Results 

We have checked the one-dimensional numerical code 
by recalculating growth rates and critical magnetic 
Reynolds numbers from the literature. In particular, we 
have considered the case of (infinitely extended) insulat- 
ing or perfectly conducting media surrounding the flow 
(for these tests we took d = 0.3 ro, and (Tsh = O.OltJfl or 
(Tsh = lOOfTfi). Using the velocity profiles corresponding 
to the models by or , we reproduced the correspond- 
ing critical magnetic Reynolds number with an accuracy 
of 5% or better. For the cases with analytical solutions 
Qil^i the accuracy of our numerical results is better than 
0.1%. 

In Ref. [9], it was demonstrated that the radial profile 
of the longitudinal velocity as measured in a water 
experiment is reasonably well approximated by 



v,ir)^U 



cosh(^) — cosh(^r/ro) 
cosh(^) — 1 



(24) 



for ^ « 18. Accordingly, we will use this parameterization 
throughout this paper, together with 



u;{r) = X 



ro 



(25) 



and treat U, ^ and x &s free parameters. One advan- 
tage of the profile (|24|l is that it provides a smooth in- 
terpolation between the laminar Poiseuille solution (for 
f — > 0) and rigid-body motion (for ^ — > oo). The latter 
limit corresponds to Ponomarcnko's model; in practice, 
for — 100, d — 5ro and (Tsh = 1, the critical Reynolds 
number differs from Ponomarcnko's solution by less than 
0.1%. 

The threshold for dynamo action depends on the con- 
ductivity and the thickness of the shell and on the veloc- 
ity profile. In Fig. Owe show the dependence of the crit- 
ical magnetic Reynolds number Rm* on the shell thick- 
ness d and the velocity profile parameter ^ for a mode 
that is close to the easiest excitable one. In the first case 
(Fig. |3^) the conductivities of fluid and shell are equal, 
in the second case (Fig. |31d) the shell conductivity is five 
times higher than that of the fluid, which approximately 
corresponds to the combination Najg / Cu. 

For the case (Tsh/cfl = 1 shown in Fig.|3^, the qualita- 
tive dependence of the dynamo threshold on the profile 
parameter ^ is different for different values of the wall 
thickness. In the limit of a thin wall {d 0), Rm, 
increases monotonously with This is explained by 
the decreasing capability of the shell to 'anchor' mag- 
netic fields lines, which can partially be compensated by 
a wider shear region in the fluid, but becomes vital as 
the flow approaches rigid-body motion for large ^. For a 
thicker shell, we find a very shallow minimum of Rm^, (^). 

For the case of a highly conducting shell, (Tsh = Scth, 
shown in Fig. ^jp, the situation is somewhat different. 
As ^ increases and thus v approaches rigid motion, the 
threshold decreases, because now anchoring of the field 
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FIG. 3: Dependence of the critical magnetic Reynolds number 
Rm» on the profile parameter ^ and the shell thickness d. The 
graphs show isolines of Rm, (d, ^) for two conductivity ratios: 
a) (Tsh = o"fl and b) ash = 5(7fl. All data are for the mode 
m = 1, k = 1/ro, which is close to the fastest growing mode. 



lines is always given (for the values d > 0.1 ro consid- 
ered here), and the shearing of field lines is maximized 
by approaching the discontinuous velocity profile. For 
fixed f > 7, the curve Rm,t((i) possesses a minimum at 
d w 0.3-0.4. Thus for the experimentally interesting val- 
ues ^ « 10-20, there is an optimal thickness d of the 
conducting shell. 

Detailed investigation of the linear dynamo problem 
in Ref. |^ has shown that in an optimal experimental 
setup the electric conductivity of the shell should indeed 
be approximately five times larger than that of the liquid 
sodium. Provided the shell is thicker than about 15-20% 
of the inner torus radius rg, the actual shell thickness 
has little effect on the excitation properties of the mag- 
netic field. For mechanical reasons, the shell should be 
kept thin (about 10% of the inner torus radius), but for 
the numerical methods applied in this paper, it is more 
convenient to consider a slightly thicker shell (typically 
about 30% of the inner radius). 

All these results were derived for CToxt = 0, i.e., an ex- 
ternal insulator surrounding the shell, which corresponds 
well to the experimental setup. In the three-dimensional 
simulations presented in Sec. IIVI however, we can only 
approximate such an insulator by setting (Text to a low 
but finite value. A discussion of the error involved by 



this approximation will be given in Sec. IIV Al 

IV. THREE-DIMENSIONAL MODELS 

If the velocity field is time-dependent, Eqs. jSJ and® 
can no longer be reduced to an eigenvalue problem, 
and we are faced with a Cauchy problem for the time- 
evolution of the magnetic field. Numerical solutions of 
this problem for the case of 2;-independent velocity fields 
have been presented in j0| and demonstrate that the 
screw dynamo should in principle work for the planned 
dynamo apparatus, which has an outer (torus) radius 
R = 0.4 m and an inner (pipe) radius tq = 0.12 m. To 
make quantitative predictions, however, the longitudinal 
dependence of the velocity field needs to be taken into ac- 
count, since only after a time comparable to the braking 
time will all of the fluid be in helical motion and thus be 
able to generate magnetic field (sec e.g. Fig. [7||. In this 
Section, we present results for this time- and z-dependent 
flow, obtained with a three-dimensional MHD code. We 
still make the approximation of (periodic) cylindrical ge- 
ometry and assume the velocity to be a given function of 
r, z and t. The problem addressed now is thus a three- 
dimensional kinematic dynamo problem with space- and 
time-dependent velocity. 

The equation we are solving is the induction equation 
in the form 

OA 

— = V X B - ry/ioMiJ + ?/oV V • A 
ot 

= (v + 7?Vln/ir)xB + ?7V2A- (r;-7?o)VV -A (26) 

for the magnetic vector potential A, from which the mag- 
netic fiux density B and the electric current density j 
are derived as B = V x A, and j = V x (B//xo/ir) — 
(V X B — V In^r X H)/ fiofir- Equation (|26() corresponds 
to the gauge in which the vector potential A and the 
scalar (electric) potential $ are related by 

77oV-A + $ = 0, (27) 

where the constant parameter rjo (introduced for purely 
numerical reasons) is arbitrary and was chosen equal to 
the magnetic diffusivity 779 of the fluid. 

As before, we use the parameterization (I24f) . H25|) for 
the radial profiles of and lu, where x is either 1 (in 
Sec. IIV A|) or z-dependent and determined from a more 
sophisticated model (Sec. IIVB|) . All dimensional results 
in this section refer to the fiducial experimental appara- 
tus ^ with the parameters given in Tabled 

For numerical reasons we have smoothed the radial 
profile of magnetic diffusivity ri{r); the resulting pro- 
file is shown in Fig. 0] the ratio 77sh/'7fl is equal to 0.2. 
We embed the cylinder in a region of enhanced mag- 
netic diffusivity; while an insulating medium corresponds 
to 77 = 00 and very large values of magnetic diffusivity 
are thus desirable, numerical requirements limit the val- 
ues of 7? strongly. We found that r?ext = Sryfl still re- 
sults in tolerable time steps, while providing already a 



6 



TABLE I: Parameters and results for the different numer- 
ical calculations presented in Sec. II VI Parameters common 
to all models are: torus radius R = 0.4 m (resulting in a 
cylinder length Lz « 2.5m), pipe radius ro = 0.12m, outer 
shell radius ri — 0.16 m, initial angular velocity of the torus 
(before braking) Qq — 310 s~^, and the magnetic diffusivities 
Vab = 0.016 m^/s, r/ext ~ 0.4m^/s. Tj, denotes the braking 
time. The amplification factors Fnct and Fmax are defined in 
Eqs. 



Label 


Vfi (mVs) 


initial field 


n (s) 


Lnct 


Tmax 


Run 1 


0.08 


random 


0.1 


87 


4.4x10^ 


Run lb 


0.08 


k = ka 


0.1 


1.2x10^ 


2.1x10'' 


Run 2 


0.04 


random 


0.1 


6.8x10"* 


1.0x10** 


Run 3 


0.08 


random 


0.2 


< 1 


1.9x10^ 


Run 4 


0.04 


random 


0.2 


860 


3.2x10" 



good approximation to the case of a surrounding insula- 
tor (cf. Sec HvH) . 

We use a numerical scheme that is of 6th order in space, 
and perform 3rd-order explicit time-stepping. Despite 
the cylindrical geometry, we use a Cartesian grid, which 
avoids the special treatment the axis would otherwise re- 
quire. The same approach was used with a similar code 
in to model nonlinear screw-dynamo action in spiral 
Couette flow. Our boundary conditions are periodic in 
the vertical direction (corresponding to the model em- 
ployed in Sec. illljl . In the horizontal direction, on the 
Cartesian boundaries of the high-diffusivity region, we 
require the magnetic field to be normal to the bound- 
aries, 

A^=0, ^A||=0, (28) 

where _L and || indicate the directions normal and parallel 
to the boundary, and d/dn denotes the normal derivative. 
Conditions H28|) imply the absence of currents across the 
boundaries which makes it a qualitative local approxima- 
tion to the case of a surrounding insulator. 

We start with a smoothed random magnetic field. By 
this we mean that the vector potential is set to a zero- 
correlated random field in all grid points within r < ro 
and to zero otherwise; we then let the magnetic field dif- 
fuse for a fraction w 0.2 of the diffusion time r^/rj^, which 
reduces the amplitude of the high-wavenumber modes 
(which decay quickly anyway) and causes the field to 
slightly extend into the shell and, very weakly, into the 
external medium. 

For the following discussion, it is helpful to characterize 
the modes by their longitudinal wave number k. Strictly 
speaking, this is only appropriate for ^-independent ve- 
locity fields like the tests given in Sec. lIV Al or towards the 
end of the dynamical calculations. In these cases, indi- 
vidual magnetic modes will have wave numbers fc„ = nki 
which are multiples of fci = 2-k/Lz. However, in many 
cases one can clearly count the number n of reversals of 




0.00 0.05 0.10 0.15 0.20 

r [m] 

FIG. 4: Radial profile of magnetic diffusivity as used in the 
numerical calculations for the resolution Sx = 0.011 m. The 
region < r < ro corresponds to the fluid; ro < r < ri 
defines the shell, and r > ri represents a poorly conducting 
medium surrounding the shell. Note that for the tests at 
5x = 0.0056 m, the profile was steeper and more step-function 
like. 



Bz along the cylinder, and we will identify this number 
n with a corresponding wave number fc„ . The azimuthal 
wave number of any relevant mode is m = ±1 (except 
for Sec lIVCll : for clarity, we will refer only to m = 



A. Steady velocities 

To test the code and to get an estimate of the error in- 
troduced by the approximations discussed above (namely 
the approximation of the surrounding insulator by a re- 
gion of enhanced magnetic diffusivity, and the smoothed 
radial 77-profilc as shown in Fig.^, we have compared re- 
sults of the three-dimensional code with those from the 
eigenvalue problem in Sec. lIIII for steady flows. 

Figure El shows the growth rate 3^7 of the magnetic 
field obtained by the three-dimensional code (diamonds 
and asterisks) as a function of Rm. The data points 
are in good agreement with the reference curve, which 
was obtained with the one-dimensional code and has an 
estimated error of < 1%. 

The validity of approximating the insulating exterior 
region by a low-conductivity medium can be explicitly as- 
sessed from Fig. El where the critical magnetic Reynolds 
number Rm* is shown as a function of (Tgh/coxt- An in- 
sulator is represented by the limit crsh/ccxt — > 00. It is 
evident that CTsh/ccxt = 25 results in a generation thresh- 
old quite close to that for dcxt 0, and the agreement is 
particularly good for ash — 5crfl. This is intuitively clear 
because the magnetic field is located around the interface 
between the liquid and the shell, with low amplitude at 
the outer shell surface. 



7 



0.06 
0.04 
0.02 
^° 0.00 
-0.02 
-0.04 
-0.06 



10 



C/[m/s] 
100 



100 



1000 




Id 

O O O 5,i:=0.011 
)^ )K yn 5a:=0.0056 



10 



100 
Rm 



1000 



FIG. 5: Kinematic growth rate K7 as a function of mag- 
netic Reynolds number for d/ro = 0.333. Solid line: one- 
dimensional model with 200 radial points. Diamonds (o) and 
asterisks (*): values obtained with the three-dimensional code 
at resolution 5x — 5y = 0.011 and Sx = 5y = 0.0056, respec- 
tively. The longitudinal wave number is fc = fcs in all cases, 
i.e., the longitudinal extent of the pipe is three wavelengths of 
the magnetic mode. The second abscissa shows the velocity 
in an apparatus characterized by ro = 0.12 m, 77 = 0.08 m^/s. 
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FIG. 6: Generation threshold Rm* as a function of the con- 
ductivity ratio (Jsh/o"oxt as obtained with the one-dimensional 
model described in Sec. IIIII The solid lines are for cTsh = cta, 
the dashed line for Ush = 5crfi . Thin lines correspond to a thin 
shell {d = 0.15 ro), thick lines to a thicker shell {d — 0.3ro). 
All values are for ^ = 10. 



B. Dynamo action in a time-dependent helical flow 

Having confirmed the accuracy of our numerical code 
for the case of constant velocities, we are now in a 
position to investigate the problem for a time- and z- 
dependent velocity field. We consider only the case of 
one single diverter in the channel, because it turns out 
that additional diverters, while shortening the transition 
time after which the flow is fully helical, have a negative 
impact on the maximum flow velocity and accelerate the 
decay of the velocity U{t). The net effect of increasing 
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FIG. 7: Evolution of the velocity structure with time. The 
braking time is Tb = 0.1, and only one diverter is used. Top: 
velocity as a function of time. The three curves represent 
the longitudinal velocity on the axis, U = Vz{r—0) (solid 
line), and the angular velocity on the axis (multiplied by ro), 
roLo{r—0), downstream and upstream of the diverter (dashed 
and dashed-dotted line, respectively). Bottom: Angular ve- 
locity uj{r—0) along the cylinder for different times. 



the number of diverters has always been found to be un- 
favorable for the dynamo. 

We take the longitudinal velocity Vz (r, t) and angu- 
lar velocity uj(r,z, t) from the hydrodynamical model de- 
scribed in Ref. [3, together with the radial dependence 
H24|l . (|25|l . At the diverter position, we slightly smooth 
the velocity field to avoid discontinuities. We first con- 
sider a short braking time of Tb = 0.1 s. Figure [7| shows 
the time dependence of Vz and uj on the axis. The maxi- 
mum angular velocity is reached only after the end of the 
braking phase and the rotating zone needs an additional 
« 0.05 s to fill the whole torus. 

Figure shows the time-dependence of the root- 
mean-square magnetic field for Run 1, compared to the 
'optimistic' and 'pessimistic' extrapolations from longi- 
tudinally uniform models: the 'pessimistic' model is ob- 
tained by adopting the angular velocity upstream of the 
diverter for any z, while the 'optimistic' model corre- 
sponds to using the downstream value 0- For < 
t < Tb = 0.1s, the initial field decays towards a sim- 
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FIG. 8; Root-mean-square magnetic flux density as a func- 
tion of time for three-dimensional models. For comparison 
with the full model (labeled 'z-dependent'), the results of two 
^-independent simulations ('optimistic' and 'pessimistic' are 
also shown. The spatial resolution here and for all following 
graphs is Sx = 5y = 0.011, 5z = 0.084. a) Run 1 (starting 
with a random fleld). b) Run lb, starting with a clean k = kg 
mode. 



pier structure (larger scales). During the time interval 
0.1s < t < 0.2 s, the resulting mode is restructured into 
TO = 1, fc = fcs, which is eventually the fastest growing 
mode. We can define the net and maximum amplification 
factors as 



maxS,.ms(0 

Srms(O) 



max B,- 
t 



min Bi.ms(t) 



r,, 



(29) 



87, 



The corresponding values for Run 1 
r„iax = 4.4 xlO^ (see Table ITJ. 

The geometrical evolution of the magnetic field struc- 
ture is illustrated by Fig.|Hl where isosurfaces of the mag- 
netic flux density |B| are shown at eight different times. 
Note how the initial mode to = 1, A; = (the slowest 
decaying mode with vanishing vertical net magnetic flux 
<I>m = J Bz dx dy in a non-helical flow) is transformed 



into k = ki and eventually k = k^. 

Given that modes with different vertical wave num- 
bers k evolve approximately independent from each other 
(since the flow always has some z-dependence, there is 
some mixing between the modes, however), it is not sur- 
prising that the field k = ki which is dominant at < = 0.15 
does not provide a good seed field for the later growth 
of the final mode k = k^. This highlights the impor- 
tance of the initial field configuration for the net growth 
of the magnetic field strength during the experiment. To 
further investigate the situation, we show in Fig. the 
growth curve from a different simulation Run lb, with 
identical parameters, but using as initial field the final 
fleld (at t — 1.5 s) from Run 1. During the flrst 0.15 s, 
the fleld decays, since the flow is not helical everywhere. 
After t = 0.15 s, the mode k = k^ grows by a factor 
of 2.1x10^ (maximum growth), and the net growth is 
1.2x10'^. Thus, the choice of the initial field can have a 
dramatic effect on the amplification factors. 



C. Enhanced magnetic permeability 

For a paramagnetic or ferromagnetic fluid, the mag- 
netic diffusivity rja is 



1 



(30) 



which shows that one possible option of increasing the 
magnetic Reynolds number is to increase the relative 
magnetic permeability fir- In an experimental setup like 
that of the Perm experiment, where there are no movable 
parts within the fluid, such an enhancement of fir can be 
safely achieved by adding ferromagnetic particles to the 
fluid. Direct measurements of the effective permeability 
of two-phase liquids indicate that /ij- can be at most « 2 
if reasonable flow properties are to be maintained [T^ . 

To assess the consequences of such an increase in Rm, 
we have carried out as Run 2 a simulation with fi^fi — 2, 
which implies a two times lower value of 779 (while r]sh 
and 77cxt are fixed to their previous values). The pro- 
file fJ.r{r) was a smoothed step profile analogous to the 
one shown in Fig. ^ but with only one step at r = tq. 
Figure 1101 shows that the net and maximum amplifica- 
tion factors are now increased by about three orders of 
magnitude. The main effect is not a faster growth of 
B, but rather a reduced decay of the initial fleld, fol- 
lowed by a much prolonged growth phase. These obser- 
vations can be understood from Fig. as follows. Since 
the dependence 5R7(Rm) is non-monotonous, doubling 
the magnetic Reynolds number will not necessarily in- 
crease the growth rate, but in fact reduce it if Rm > 40. 
On the other hand, the critical flow speed (correspond- 
ing to Rm « 21) is two times lower than for fi^fi — 1, 
and hence the flow is supercritical for a much longer time 
span. 

The fourth line in Fig. ^| shows an artificial z- 
dependent run, where the 'paramagnetic pumping' ve- 



FIG. 9: Structure of the magnetic field for different times of Run 1. From left to right, the times are 0.0, 0.05, 0.1, 0.15, 0.2, 
0.3, 0.4 and 0.5 s; the braking time is Tb = 0.1 s. The surfaces are isosurfaces of the magnetic field strength (red: Bz > 0, blue: 
Bz < 0), the lines are magnetic field lines. The diverter is located at the bottom and the flow is directed upwards. 
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FIG. 10: Like Fig.|H| but for Run 2, i.e. for a two times lower 
magnetic diffusivity (achieved by enhancing the magnetic per- 
meability to firei = 2). Note the different time interval plotted 
and accordingly the much longer growth phase compared to 

Run 1. The fourth line ( ■ •) shows the results for a model 

where the 'paramagnetic pumping' velocity Vp was artificially 
set to zero. 

locity Vp was set to zero; this would correspond to the 
case where rja is reduced by enhancing aa, rather than 
/ir,fl. The comparison shows that the pumping term is 
indeed important for the field growth, since without it 
the net growth would be reduced to about 3400. 

Another interesting finding is that for Run 2 the rms 
magnetic field for the full model is no longer contained 
in the interval spanned by the 'optimistic' and the 'pes- 
simistic' variants. Rather, the rms field for the 'pes- 
simistic' run overtakes the z-dependent one and closely 
approaches even the 'optimistic' run. This unexpected 



behavior is connected to the presence of several grow- 
ing modes and can be understood by a closer look at 
the modal structure of the solutions. In Fig. we show 
the time evolution of individual modes identified by their 
vertical wave number k. Strictly speaking, only for z- 
independent velocity profiles (i.e. for the 'optimistic' and 
'pessimistic' models) each dynamo mode will be char- 
acterized by a unique value of k (together with the az- 
imuthal wave number m and a radial one), but spectral 
analysis for the z direction is a very helpful tool even if 
this condition is not satisfied. 

As can be seen from Fig. ^2 in the 'optimistic' model 
the mode k — k^ (and m = 1) dominates, while for the 
full and 'pessimistic' models k — kj^ (and m = 1) is the 
dominating mode. For most of the time, the growth rate 
of the modes and ^4 are comparable, which indicates 
that a wave number between k^ and fc4 would be optimal, 
but is excluded by the geometrical setup. Other growing 
modes include k — kj, m — 2, which was never encoun- 
tered for Run 1, where the magnetic Reynolds number 
was two times lower. 

Since the magnetic diffusion time rl/r] = 0.36 s is at 
least comparable to the evolution time of the flow, the 
growth of the modes is never just determined by the 
growth rate for the current value of Rm, but rather in- 
volves the history of the given mode. For the full, z- 
dependent model, an additional important factor is the 
action of the diverter, which for all times (but particu- 
larly during the first 0.2 s) introduces a z-dependence of 
the velocity field and thus mixes energy from the domi- 
nating mode k — kji into other modes. This is the reason 
why in Fig. Illb the different modes show quite similar 
behavior after t k, 0.3 s. Although the same effect will 
occur for /^r = 1, too, the resulting energy loss from 
the dominating mode will be weaker there, because the 
shorter magnetic diffusion time allows the leading mode 
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FIG. 11: Evolution of Fourier modes for Run 2. a) 'optimistic' 
model; b) full model; c) 'pessimistic' model. Shown is the 
energy Ek,ra of the modes normalized such that the value 
at f = is 1. The oscillations of the mode fc = fca in the 
second plot are connected to a change in the radial structure, 
i.e. most probably due to different radial modes. 
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FIG. 12: Like Fig. |H] but for different parameters, a) Like 
Run 1, but with a longer braking time Tb = 0.2 (Run 3). b) 
Like Run 2, but with a longer braking time Tb = 0.2 (Run 4). 



D. Slower braking of the torus 

To assess the effect of a longer braking phase, we have 
carried out calculations with Tb = 0.2 s; the results are 
shown in Fig. 112b for [i^ ~ \ (Run 3) and in Fig. 112b 
for /if = 2 (Run 4) . In the case of liquid sodium without 
admixed ferromagnetic particles, the maximum growth 
is diminished to about 100 and the net growth is less 
than 1, i.e., at the end of the experiment the magnetic 
energy is lower than it was for the seed field. Only with 
the enhanced value /ij- = 2 wc obtain net growth, which 
is now larger than in Run 1, but still significantly lower 
than in Run 2. 



to better adjust to the z-dependence of the velocity field. 

As a result of these effects, drawing conclusions for the 
full problem from the simple 'pessimistic' and 'optimistic' 
models can be quite problematic. 



V. CONCLUSION 

The results presented here confirm earlier estimates ac- 
cording to which the planned Perm dynamo experiment 
is realistic and will be able to yield field amplification 
factors of about 10'^ or more. A thin, highly conducting 
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shell is crucial for the dynamo process and its role is well 
understood. 

Short braking times are necessary for the dynamo, and 
a time of Tb = 0.1 s as is intended for the Perm experi- 
ment will be sufficient. 

Enhancing the magnetic permeability by adding fer- 
romagnetic particles to the liquid sodium would further 
enhance the amplification factor, but this is not crucial 
for the success of the experiment. For a longer brak- 
ing time Tb = 0.2 s, however, enhancing /Xj- is required 
to obtain net growth of the field at all. The enhanced 
magnetic Reynolds numbers for /ij. = 2 cause a number 
of modes to grow (which complicates the analysis) , and 
enhances the interaction of different modes due to the 
inhomogeneity introduced by the diverter. 

The final amplification factor strongly depends on the 
initial magnetic field configuration, and choosing a suit- 
able seed field can be vital for obtaining good results. 
For example, relying on the terrestrial background field 
may be a bad choice, as a uniform field penetrating a 
torus has only components m = 0, k = ±fci and m = 1, 



k = while a good seed field should have a significant 
amount of energy in the optimal mode m — I, k — k^. A 
sophisticated arrangement of small permanent magnets 
may be able to provide such a field, but one should even 
consider an arrangement of coils that allows for a net cur- 
rent through the inner part of the torus. This question 
certainly requires more detailed investigations. 
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